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Abstract 

We apply the Nozieres-Gallet dynamic renormalization group (RG) scheme to a continuum 
equilibrium model of a d-dimensional surface relaxing by linear surface tension and linear surface 
diffusion, and which is subject to a lattice potential favoring discrete values of the height variable. 
The model thus interpolates between the overdamped sine-Gordon model and a related continuum 
model of crystalline tensionless surfaces. The RG flow predicts the existence of an equilibrium 
roughening transition only for d = 2 dimensional surfaces, between a flat low-temperature phase 
and a rough high-temperature phase in the Edwards- Wilkinson (EW) universality class. The 
surface is always in the flat phase for any other substrate dimensions d > 2. For any value 
of d, the linear surface diffusion mechanism is an irrelevant perturbation of the linear surface 
tension mechanism, but may induce long crossovers within which the scaling properties of the 
linear molecular-beam epitaxy equation are observed, thus increasing the value of the sine-Gordon 
roughening temperature. This phenomenon originates in the non-linear lattice potential, and is 
seen to occur even in the absence of a bare surface tension term. An important consequence of 
this is that a crystalline tensionless surface is asymptotically described at high temperatures by 
the EW universality class. 
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I. INTRODUCTION 



The dynamics of growing surfaces [0, || has attracted great interest during the last decade. 
This is due both to the practical implications for the control of film quality in thin film pro- 
duction techniques, and to the fundamental questions it raises in areas of Physics such as 
spatially extended systems in the presence of fluctuations 0, or scale invariance in non- 
equilibrium systems It has been observed that a surface growing in the presence of 
fluctuations quite generically exhibits time and spatial scale invariance properties. Specifi- 
cally, if h{r,t) denotes the surface height at time t above a (i-dimensional substrate position 
r G M.'^, the height- difference correlation function G{r,t) = {{h{ro + r, t) — /i(ro,t))^) grows 
as G{r,t) ~ ^2^/2 for t^/^ r and scales as G{r,t) ~ r^" for t^^^ 3> r, where a and z are 
referred to as the roughness and dynamic exponents, respectively. These scale invariance 
properties imply that such rough surfaces share many properties with dynamic critical phe- 
nomena, which allows one to obtain useful information on the former from studies on the 
latter, and vice versa. A particularly interesting example is provided by the overdamped 
sine Gordon (sG) model, which describes the equilibrium fluctuations of a crystal surface 
and the features of its roughening transition [^: 



fi~'— = F + uAh sin + v/2T/i-ir7 r,t . 1 

ot ct± \ a± J 



In Eq. (0), A is the substrate Laplacian, fi is the surface mobility, F, u, V and are positive 
constants, T is temperature, and rj is a zero mean, Gaussian white noise with correlations 
{ri{r,t)ri{r',t')) = 6{r — r')6{t — t'). The surface morphology of a crystal surface described 
by the sG model thus results from an interplay between linear surface tension, described 
by the term with coefficient in (|1|), a periodic potential favoring discrete interface values 
h = n a± with n G Z, and thermal fluctuations. Moreover, using the same arguments as 
in Q, for a non-zero value of the homogeneous driving (flux of aggregating particles) F in 
(H), an alternative interpretation of the sG surface is as one minimizing surface area (for 
small values of the surface slope), subject to a flux of aggregating particles such that growth 
events of the height take place in integer values. In equilibrium {F = 0), the sG model is 
well known § to feature a roughening transition when d = 2: for temperatures 
above a critical value T > T|P, the lattice potential is irrelevant and the surface is rough, 
in the sense quoted above for the height-difference correlation function. The corresponding 



exponent values are those of the Edwards- Wilkinson (EW) equation |]TT|, which is simply 
the 1/ = limit of Eq. (1), namely 



2-d 

"EW — — - — , — 2. (2) 

In the case d = 2 considered, (||) amounts to oew = 0, implying G{r,t) ~ logr for r ^ 
t^/^. For temperatures smaller than T|P, the lattice potential dominates the large scale 
properties of the surface setting a finite correlation length ^, beyond which the height- 
difference correlation function G{r > ^,t) attains a constant value and hence the surface 
is flat. The sG roughening transition is of the Kosterlitz-Thouless class, and hence the 
continuum sG model is related to important discrete models such as the discrete Gaussian 
and the models P, |T3l . 



As mentioned above, a practical domain which presents many instances of growing sur- 
faces is the area of thin film production by techniques such as e.g. Molecular-Beam Epitaxy 
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(MBE). In many MBE conditions |jT^, the main relaxation mechanism on the surface is sur- 
face diffusion, rather than surface tension. In this case, and again considering the additional 
effect of a periodic lattice potential, one is led naturally to the following model of a growing 
surface henceforth referred to as the xMBE model 



2W (2Txh 
sin 



v/2T/i-ir/(r,t). 



(3) 



In Eq. (0), K is a positive constant, in which case the corresponding term in the contin- 
uum equation does provide — again for small surface slopes — a surface diffusion relaxation 
All the other terms in Eq. have the same meaning as in (|I|). Specif- 

Thus, for F ^ 0, Eq. (g) 
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mechanism 

ically, F is related to the flow of adatoms onto the substrate, 
provides a non-equilibrium description of a growing surface which roughens under the effect 



the surface described by Eq. 



of thermal fluctuations. In particular, as shown in 
initially displays RHEED (Reflection High Energy Electron Diffraction) oscillations, akin to 
those observed experimentally [|r^. For F = 0, describes the equilibrium fluctuations at 
temperature T (we will consider a unit Boltzmann's constant) of a surface minimizing the 
energy functional (Hamiltonian) 



d'^r^|(A/.)2 + l^ 



1 — cos 



l-nh 



(4) 



Therefore, such a surface can be also thought of as minimizing surface curvature to linear 
approximation, thus having in principle zero surface tension, and to favor values of the height 
which are integer multiples of a^. Note in this respect that (^ does not have the form of a 
continuity equation for the surface height, due to the form of the nonlinear and noise terms 
p!8| . Hence, from the point of view of applications to MBE growth, the continuum xMBE 
model (^ might be relevant for those situations in which nonconserved noise is expected to 
play a role, e.g. for length-scales larger than the typical diffusion length [|TU]. On the other 
hand, in a similar way as the sG equation can be seen as a continuum description of the 
discrete Gaussian model, Eq. (^) is a natural candidate for the continuum description of the 
discrete Laplacian roughening (Lr) model pOf on the square lattice |^T|]. The Lr model has 
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been used to describe two-dimensional defect melting and as a model of tensionless surfaces 
such as membranes [|^, Actually, Eq. (|) can indeed be obtained 
limit of the Lr model, albeit in a non-unique fashion. 

151 , p5[] show an equilibrium roughening tran- 



Numerical simulations of Eq. (^ for d = 2 



sition, similar to that in the sG model. In the case of Eq. (|^), the exponents characterizing 
the rough high temperature phase are numerically consistent with those obtained for the 
V = Q limit of Eq. (|) —the so called linear MBE (/MBE) equation p6|, namely 



A-d 



«iMBE 



^ZMBE 



(5) 



which in d = 2 implies azMBE = 1- However, given that in the sG system the periodic poten- 
tial is known to contribute a correction to the surface tension term upon renormalization, 
the same phenomenon is expected to occur in the xMBE model, in which the bare surface 
tension is zero. Since surface diffusion is irrelevant in the presence of surface tension [0, ^ , 
this argument together with the available numerical evidence |1T5|, EH] led to the prediction 



that the scaling properties of the xMBE model in the high temperature phase are the same 
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as those of the sG modeL To date, and although more recent numerical data confirm 
the existence of crossover behavior in the critical properties of model (^), direct numerical 
verification of EW scaling properties seems hard to achieve. Moreover, previous analytical 
studies , which employed a variational mean- field analysis successfully applied in for 
the analysis of the sG roughening transition, have failed to reproduce EW scaling at high 



temperatures for the xMBE model. Specifically, the results obtained in [E3] include an upper 



critical dimension dc = 4, below which model (0)-(^) displays a phase transition between 
a low temperature fiat phase and a high temperature rough phase whose scaling behavior 
is that of the /MBE equation. However, the transition taking place is of first order for all 
substrate dimensions, d < dc, while the numerical simulations of [1T5| , p5[] strongly suggest 



that, at least for d = 2, the transition is of a continuous type. Besides, the generation of a 
surface tension term by the lattice potential referred to above is due to non-trivial coupling 
between different surface Fourier modes, which is unsufficiently accounted for by the vari- 
ational mean-field description employed in [^. Therefore it is natural to try improve on 
the mean-field approximation of Ref. ||2^ and explore the scale invariant properties of the 
system near the phase transition point through the use of the renormalization group (RG). 

In this paper we consider the following generalized continuum model of surfaces subject 
to a periodic lattice potential: 

IJ,^\- = F + uAh - nA^h sm ' 



dt a± \ a± 

+ y/2Tfi~^r]{r,t). (6) 

This equation obviously features both Eqs. (0) [sG model] and [xMBE model] as special 
cases. Moreover, the linear limit V = ||T^ ^ of this equation has been observed to 
accurately describe growth experiments of copper aggregates by electrochemical deposition 
in the presence of organic additives For the sake of simplicity, in this work we will 



consider this generalized model in the absence of driving {F = 0), in which case Eq. 
describes the equilibrium fiuctuations of a surface with Hamiltonian 



E,{u,K,V) = Jd'rS^'^iVhf+'^iAhr + V 



cos 



27Th 

a± 



(7) 



In order to study the critical properties of the xMBE model in d substrate dimensions, we 
will extend the dynamic RG approach devised by Nozieres and Gallet (NG) for the d = 2 sG 
model [§] (see also a detailed account in 0) to the generalized system (|]). There are two 
reasons for our pursuing this approach: (i) as anticipated above, in the RG analysis of (H)-® 
a finite (non-zero) surface tension term needs to be allowed for, given that it is generated in 
any perturbative scheme even if its bare amplitude is zero; (ii) a static renormalization group 
study of the equilibrium system (^ is ill-defined in some parameter ranges due to divergent 
integrals ||2^, similarly to the sG case p. Still, the static RG study will provide us, via the 
appropriate generalization, with the correct expansion of the model non-linearity in terms 
of relevant operators through the use of Kadanoff's operator product expansion (OPE) []32| , 
as was accomplished in ^ ^ for the sG model. In any case, the results to be obtained 
from the dynamic RG study that follows will also cover the case of a system minimizing 
both surface area and surface curvature, and will in particular allow us to analyze how the 
standard sG roughening transition is modified when an additional surface diffusion term 



4 



is considered. We will finally consider the renormalization properties of the xMBE model, 
which corresponds to a specific choice of bare parameters within this generalized framework, 
and will compare the conclusions obtained with those from both the variational approach 

2^ and numerical simulations [^, |2^. To our knowledge, ours is the first RG approach 
to the xMBE model as formulated by ®-(^, and it may contribute to the elucidation 
of the existence and nature of the phase transition in this and related systems, as the Lr 
model. There also exist static and dynamic RG studies of similar systems. Specifically, the 
equilibrium properties of a model which is different from (|^)-(^) but is believed to provide 
the continuum description of the Lr model on the triangular lattice, have been analyzed in 

3^ , and its dynamical properties have been obtained in and references therein. Within 



the rough MBE surfaces context, the dynamic properties of the conserved sG model have 
been studied both under conserved [111, E 



T8| , p% noise. Finally, a 



I , pql and non-conserved 
Langevin equation beheved to describe a restricted curvature model has been analyzed 
in by using RG techniques. 

This paper is organized as follows. In Sec. || we apply the dynamic RG scheme of NG 
to Eq. (^. The parameter flow thus obtained for this generalized model is studied in Sec. 
Ill A| as a function of the substrate dimension d. The special limit of Eq. corresponding 
to crystalline tensionless surfaces, Eq. (|^), requires additional considerations of a technical 
nature, and is deferred to Sec. |III B| . Finally, Sec. |I^ is devoted to further discussion of 
the results obtained in the previous sections, and to summarizing our conclusions. We 
additionally provide two appendices. In Appendix ^ we detail, following ||, the OPE 
which is needed in the dynamic RG in order to perform the appropriate expansion of the 
lattice potential in (^) into relevant operators. Appendix ^ closes with a discussion of the 
specific way in which the roughening transition of the sG model generalizes into that to be 
obtained in Sec. 00 for 



II. DYNAMIC RENORMALIZATION GROUP ANALYSIS 

This section is devoted to the analysis of Eq. in the equilibrium case F = 0, em- 
ploying the dynamic RG scheme of NG @, 0. In this scheme, a coarse graining procedure 
is performed over the microscopic modes of the noise term. Namely, the noise is split into 
two statistically independent parts, rj = fj + Si], such that the total noise power spectrum 
is the sum of the corresponding contributions. Here 6r]{r) = d^k e^'*^'"?7(k), where A is 
a momentum cut-off related e.g. to atomic positions on the substrate, A = e~^A with e a 
small parameter, and i7(k) is the spatial Fourier transform of the noise. Then, microscopic 
fluctuations are integrated out by defining h = {h{fj + Sri))srf and Sh = h — h, and by seeking 
an equation of motion for the thus defined long distance modes h. The result will be an 
equation with the same shape as (|^), but with new (renormalized) coefficients, which are 
sensitive to the microscopic fluctuations 6t] by the action of the nonlinearities. Specifically, 
the dynamic equations for the h and 6h modes read 

dh - - - 

/x"^— = uAh- nA'^h- {^{h,5h))sr^ 

+V2Dr]ir,t), (8) 

dt' ^ 



+V2D6r]{r,t), (9) 
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where we have defined D = Tfi~^, and we have introduced Sh) = {2TTV/a±) sin[27r(/i + 
6h)/aj_]. In order for (H) to be a closed equation in h, we need to solve for 6h in (||) and 
introduce the result into fBI). The formal solution of fOf) reads 



6h{r,t) = / dV /" dt'xo{r~r',t-t') 

J J —CO 



(10) 



where the primed notation denotes dependence on the r', t' variables, and the d-dimensional 
free propagator reads 



Xo(r,t) 



(11) 



Due to the non-linear lattice potential, an explicit solution of (|^) can only be obtained by 
performing a perturbative expansion in powers of V. Thus, defining 



we obtain 



(5/i(°)(r, t) = f dV f dt' xo(r ~r',t- t') 

J J ~oo 

xV2D67]{r',t'), 
6h^'\T,t) = / dV r dt' Xo{r-r',t-t') 



(12) 



x5h'^^^ cos 



27r/i' 



(13) 



(14) 



Using sin(a + 6) = sin a cos 6 + cos a sin 6, and within our perturbation expansion, we can 
now evaluate in (B) 



{^{K5h)) 



St) 



1W 



l-^((5/.(°))% 



a I 



47rV 



X sm 



which has the form 



with 



(27r)2-% T A'^-^eT 



(15) 



(16) 



$«(/i) 



27r 
xsin 



2ai(i/ + fi:A2) 

27r/i 

5r, sin 



27r/i 



(17) 
(18) 
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where Sd = 27r'^/^/r((i/2) is the surface area of the unit hypersphere in d dimensions, r(-) 
being Euler's Gamma function. The shape of ([TTD aheady reflects the fact that we have 
considered an infinitesimal shell of microscopic noise modes of width e. We can write (|18D 
more explicitly as 



32tt' 



I dV f xo{r-r',t-t') 

J J — OO 



(19) 



x{5h^^^6h'^'^)s^ si 



sm 



2'Kh 
a± 



cos 



2nh' 



Neglecting higher order harmonics of the lattice potential, sin (27ih/a±j cos{2'7ih' /a±) ~ 
^sm[2'7i(h — h')/a±]. Further, we can use the results for the OPE of the lattice potential 
obtained in Appendix |^, and a Taylor expansion to obtain 

327r« 



) /»oc poo 

-Sj p'-Mp / drxo(p,r) 
Jo Jo 



(20) 



. dh 

^ ~dt^2d ~ 8d{d + 2) 



2'kT 



a_L 



0(p, r, z/, K, p) 



x(5/i(°)(5/i'^°^)5^exp 

where we have defined p = r — r' and t = t — t', the ci-dimensional free propagator reads 

P 



Xo(p,r) 
and we have introduced 



(27r)'^/2 



A i.d/2 

dA;4^J,/2-i(pA:)e-('^^'^+-^^) 



P' 



d/2- 



(21) 



(27r)^-i 



A ud-l 



1 - 



T{d/2)Jd/2^i{pk) 
(pA;/2)'^/2-i 



(22) 



with Jn{x) being the n-th order Bessel function of the first kind. Finally, using the results 
of Eqs. (|l3) and (^), Eq. (0) has the form 



{<^)5hCh) = {V + 6SV)sm 



2'Kh 

2; 



+eSuAh + eSKA^h + 0{V 



(23) 



where the corrections Su, 6k,, etc. are implicitly defined. Inserting the result of Eq. (^) into 
Eq. (H), we obtain that, to order, the long distance modes h obey an equation with the 
same shape as Eq. (10), namely, 



p — = u Ah — kA n sm 

ot a± 



27rh 



(24) 
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but with new coefficients u = u + e Su, k, = n + e Sk, fi ^ = ^ ^ + e5^ ^, and V = V + e 5V. 
In order to recover the original Fourier mode cut-ofF A, we now rescale variables as 



r 

h 
t 



where 6 = e^, and we thus get 



^ dt 



r' = r/b 

h' = hb-" 
t' = b-H 



■sin I ) + V2Df]{r,t), 



a± \ a± 
kb'-\ V = 1/6^-2°, 



(25) 
(26) 
(27) 



(28) 



'z—2a—d 



with coefficients z/ = R = kb""' , V = V b^~ , ft = fi~ , and D = Db^ 

Finally, for an infinitesimal e = dl and expanding iy{s), etc., to first order in e, we obtain 
the dynamic RG fiow to V'^ order: 



du 
dl 



dn 
'dl 



dV 
'df 
dfj-^ 

~dr 

dD 

'df 

da I 



V2 5(2'0)rzy, f,) 



^2 5(4.0)/ X 



dl 



{z - 2)iy 
(27r)6-'^5rfTA'^-6 

(z - 4)k 

(27r)6-^5rfTA'^-8 

^ ' 2a^(z/ + KA2) 

(2vr)^-^5.TA'^-^ _ ,0,, 

{z-2a- d)D, 
—aa±, 



(29) 
(30) 

(31) 

(32) 

(33) 
(34) 



with 



G{p, r, z/, k) 



oo roc 



dpp 

JO 

1 _ _ 

d/c k"^/^ Jd/2-iipk) exp 





dr e ^ r"* J^/a-i (p) ^(p, r, z/, k) exp 



— T- 



■'^^(f){p,T, Z/, K,/^5) 

(2 I 



(36) 



where p = pA. It is worth noting that the fiow of the mobility p is enslaved to that of all 
other system parameters, which allows us to neglect its evolution under (PUD-(P1D in the rest 
of the paper. Moreover, in order to preserve fiuctuation-dissipation by this RG procedure 
P, 1^ we need to impose the exponent relation z = 2a+d in the fiow (p9|)- (|3^) . Note that this 



relation (termed hyperscaling |40] in the studies of kinetic roughening) holds exactly at two 
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fixed lines of (p9|)-([3^), wliicli are tlie two linear limits of namely the EW [V = k = 0, 
see Eq. (j^)] and /MBE [V = u = 0, see Eq. (|])] equations. Both systems can be interpreted 
as describing equilibrium fluctuations, governed by the corresponding Hamiltonians. We 
also remark that, as anticipated in Sec. I, it is clear from Eq. ( pQ]) that, if we consider the 
RG flow for the xMBE model (0)-(^ in which there is no bare surface tension, the lattice 
potential does generate it under RG iterarion. 

Finally, we find it more convenient to express the flow in the dimensionless variables 



Thus, we have 



X ; 

dx 
d/" 



2al 
Iff 



y = 



AW 



K 



{d - 2)x - AK 

M^\x,K) - 



+ 



Ay' 



dy 
dl 
dK 



d 

'd 



A{d + 2) 



B^^\x,K) 



(2n 



\d-l 



{d-A)K- 



2d{d + 2) X' 



B('\x,K), 



(37) 
(38) 

(39) 
(40) 



where 

B'^''\x,K) 
G{p, T,x,K) 
0(P, r,x,K) 



dp p 



(2vr)^-i Jo 



dre ''Jd/2-iip)G{p,T,x,K)e' 
dk~k'/' Jd/2-iCkp) 6-2-1(^/2-^)^^+^^*]/-, 



-24,ip,r,x,K) 



S,A 



d-2 /•! 



dk 



(27r)^-i Jo (x/2 - K)k' + Kh 



^ _ T{dl2)Jd/2-i{kp) 
(^p/2)^/2-i 



(41) 
(42) 

T[{x/2-K)k'^+Kk'']/^W^ 



Note that the convergence properties of the integral 0(p, r, x, K) defined above depend on 
the substrate dimensionality d and on the parameter values. Specifically, for d = 2 and 
the condition associated with the xMBE model (^, namely, x = 2K, the integral diverges 
logarithmically at the lower limit. This divergence originates in the fact that correlation 
functions do not have a well-defined thermodynamic limit for the /MBE model in d = 2 fLA . 
The dynamic RG flow (|38|)-(pOD just obtained generalizes that in p, ^ for the case of 



a finite surface diffusion term, and for any value of the substrate dimension d. The sine- 
Gordon case is retrieved simply by setting n = and neglecting the higher order contribution 
B'^'^^x, 0), related with the RG flow of the surface diffusion term. Actually, this integral will 
allow us to study in the next section the effect of such relaxation mechanism on the sG 
roughening transition. 



III. DISCUSSION OF THE RG FLOW 

The fixed point structure of the RG flow (^8])-(0) depends on substrate dimensionality 
d, and thus so do the critical properties of model (D-(0). In what follows we will restrict 
ourselves to dimensions d > 2 in which phase transitions are expected to occur. 
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For general values of d, there are no proper non-trivial fixed points of (p8D-(^), but 
rather regions in the (x, y, K) phase space which are invariant under the RG flow. These 
correspond to the two significant linear limits of the generalized model, namely the EW line, 
FeW) and the /MBE line, F^mbe, defined by 

Tew = {{x.y.K) -.y = K = Q,x ^Q}, 
T/MBE = {{x,y,K) : X = 2K, y = 0}. 

However, for the special dimensions d = 2 and d = A, respectively, these regions become 
lines of fixed points. 



A. Generalized model 

We start by analizing generic features of the RG flow as a function of substrate dimension 
d, in the case in which the bare values of surface diffusion and surface tension are both non- 
zero. The special initial condition (for the RG flow) in which there is no bare surface tension 
(xMBE model) will be considered in detail in Sec. |111 B . 



1. d = 2 case 



Obviously, the case which is most interesting from the physical point of view is that of 
a two-dimensional substrate. As stated above, now the invariant EW region Few actually 
becomes the well-known sG line of fixed points p. An important point along this line is 
a; = 1, where the flow of the lattice potential y changes stability, see Appendix |^. If the bare 
value of the surface diffusion k{1 = 0) = 0, as in the sG model, numerical integration of (|5BD- 
(PH) shows that the flow essentially remains on the (x, y) plane and thus completely reduces 
to that of the sG system, as seen in panel (a) on Fig. |l|: for high T the flow is towards the 
y = axis, the interface being rough and characterized by the scaling properties of the EW 
equation. For low T values, the lattice potential y grows indefinitely upon iteration of the 
RG flow, drawing the system into the T = limit for which the interface is fiat. Thus, there 
exists a roughening transition between these fiat and rough phases. Numerically, we estimate 
the critical temperature as the value characterizing the trajectory that separates between 
flow onto the y = line and flow towards increasing y. We thus obtain ~ 0.725 ± 0.05 
for V = l,a± = l,z/ = 1,K = 0, which is close to the value for the flow equations of the 
pure sG model p, |[, = 2/7r ~ 0.64. On the other hand, as we see in panels (b) and 
(c) on Fig. 0, in general a small initial surface diffusion shifts the roughening temperature 



to higher values. The reason for this effect is that, as studied in when V = 0, 

even though surface diffusion is irrelevant relative to surface tension in the hydrodynamic 
limit, an initial value k{0) ^ introduces a crossover length scale, L^il) = [fi:(Z)/z/(Z)]^/^, 
below which surface diffusion is the relevant relaxation mechanism. Once the coarse-graining 
procedure has overcome this scale, renormalizes to zero (see Appendix |B|) and the flow 
takes place on the sG (x, y) plane with a u coefficient which differs from its bare value. 
From there on the flow is effectively as that of the pure sG model. The results shown on 
panel (d) of Fig. |^ correspond to a numerical simulation of the sG model [Eq. (^ with 
K = 0,F = 0] for fi = i' = V = a± = l, and of Eq. (^ for k, = 0.5, for both of which we have 
computed the specific heat [defined as xe = {{Eg) — {Eg)"^) / {TL"^), with Eg as in (|^)] for 
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several substrate lateral sizes L and periodic boundary conditions. Clearly, the roughening 
transition temperature — which is preceded by a peak in xe as in the pure sG model [J2 — 
shifts to higher values for n ^ 0. Integrating numerically the RG flow (|38D-(pOD, we obtain 
Tf{K = 0.5) ~ 1.35 ±0.05. 



2. d > 2 cases 



In order to inquire about the upper critical dimension, dc, of model we consider 

values of d above the physical two-dimensional case. An important consequence of this is 
that the points on the EW line Few are not fixed any longer under the RG iteration. This 
already signals the final result that, in these dimensions, there is no proper phase transition. 
Rather, the only existing phase is the fiat low temperature one. Before justifying this result 
in detail, let us note that the EW equation already predicts a flat surface for d > 2 p|, |^, p!2 



and therefore when adding a lattice potential to the surface tension term in the sG equation, 
the result that = 2 is the corresponding upper critical dimension ||3^, ^ does not come 
as a surprise. In the generalized model (|§)-(0), similar results to those in the previous 
section indicate that the same value dc = 2 actually occurs: for 2 < < 4, it is clear from 
Eq. ( |40D that K{1) still decreases exponentially under the RG iterarion, and again the flow 
becomes essentially that of the sG model in the corresponding dimension, thus predicting a 
fiat morphology. In Fig. ^(a) we show as an illustration results of a numerical integration of 
(^)- (|40|) for d = 3. As is clear from the figure [for this, it is useful to note the projection of 
the RG flow lines onto the {x,y) plane], the flow is eventually towards large y values for all 
temperatures. Note also that the K{1) value may in some cases become negative, signalling 
a non-physical instability. We attribute this effect [also apparent on panels (b) and (c) of 
Fig. 3 to limitations of our approximation to the RG flow. 

The value d = 4 is marginal for the surface diffusion term in Eq. (^). In an analogous 
way to the role of d = 2 for the EW equation [see Eq. (0)], the /MBE equation predicts a 
roughness exponent a = for = 4 [logarithmic behavior of G{r,t), see Eq. (^] and a flat 
morphology ior d > 4 |1|, 0. Even though the decay of K{1) under the RG flow (|38|)-(^0[) 
might not be so fast as in smaller dimensions, it nevertheless occurs, see Fig. 0(b), with 
the result of an effective sG behavior in ci = 4 dimensions, again corresponding to a flat 
morphology. Finally, as illustrated in Fig. ^(c) for ci = 5, if ci > 4 all coefficients in the 
RG flow grow under iteration. In particular the lattice potential y increases indefinitely and 
thus the surface morphology is dictated by the behavior of the T — > limit, namely again 
there is no phase transition and the surface is flat. 



B. Crystalline tensionless surfaces 

As remarked above, within the RG flow of the generalized model (^)-(|3)) the case cor- 
responding to the xMBE model (§)-(§) simply amounts to a specific condition on the bare 
parameters, namely they lie on the plane a;(0) = 2K{0), y{0) ^ 0, thus implying z/(0) = 0. 
In this section we thus study the RG flow for the xMBE model separately for different values 
of the substrate dimension. 
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1. d = 2 case 



In the physical d = 2 case, and due to the divergencies in the correlation functions 
mentioned above, this very condition induces a trivial RG flow, since the S^"^ integrals in 
(IT) become identically zero. The linear combination x — 2K — which is proportional to 
the surface tension — , is a constant under the RG flow. Equations (p8|)-(P0|) can be exactly 
solved to show that any initial condition on the x = 2K plane tends to the origin for I ^ oo. 
For large enough but finite /, the surface scales as the /MBE equation independently of the 
value of T, and thus there is no temperature driven phase transition, which contradicts the 
results of e.g. numerical simulations of Eq. (|^) for two-dimensional substrates Thus, 



some kind of integral regularization is needed when d = 2. Here we introduce a lower 
momentum cut-off 1/(AL) in the integrals (021), (^5D, with L being a measure of the lateral 
dimension of the two-dimensional substrate. Proceeding in this way, it can be seen that 
the integrals ( ^T] ) now vanish as a power law of 1/L. In this regularization scheme, the RG 
iteration has to be stopped once e' = L. Another possible procedure like that of dimensional 
regularization can be employed with the same conclusions as those that follow, similar to 
the RG analysis of tensionless membranes p4| . 

In any case, the regularized integrals B^^\x,x/2) are no longer identically zero. This 
leads us to expect the RG flow to escape from the x = 2K plane, which would mean a 
finite surface tension has been generated under renormalization, and thus the occurence of 
temperature dependent behavior. Nevertheless, the flow may take many iterations before it 
appreciably deviates from the x = 2K plane, given that, in its neighborhood, the integrals 
i?*^") can be rather small numbers for large L values. 

In Fig. 1^ we show the numerical integration of the RG flow (|38|)-(^0D for various initial 
conditions on the xMBE plane and L = 128. Indeed, three different types of behavior can 
be distinguished. For high enough values of T, the variables x, y, K are very small and the 
RG flow escapes from the x = 2K plane very slowly. Under these conditions, the system 
flows towards the origin while featuring the scale invariant behavior of the /MBE line F/mbe- 
An example is trajectory (1) on Fig. |^. For intermediate temperatures, v ~ x — 2K becomes 
non- negligible and the flow is attracted by the (x, 0, 0) segment with x < 1, where the 
behavior is described by the EW equation, see trajectory (2) on the figure. Finally, for high 
values of T, the flow falls rapidly onto the (x, y) plane with x > 1, see trajectory (3) on the 
figure. This behavior is described by the low temperature massive phase of the sG model, 
where the lattice potential is dominant and the surface is flat, see Appendix ^. 

In order to examine this behavior more closely, we can focus on the crossover length Lx 
introduced in Sec. [Ill A[ By coarse-graining the system with the RG transformation for scales 
up to the system size {I* = InL), and if the lattice potential turns out to become irrelevant, 
we can decide whether the system scaling is of the /MBE or EW type by evaluating Lx (/*) to 
be larger or smaller than L{1*) = 1, respectively. Note that, for a zero lattice potential, the 
RG flow of Lx is easily computed from (p9| ) -(|30| ) to be given by Lx (/) = Lx (0) exp(— /), hence 
an exponential behavior of Lx(/) will be a sign of the irrelevance of the lattice potential at 
scale /. In Fig. ^ we show the RG flow of Lx (/) for the same temperature conditions as in Fig. 
^ and L(/ = 0) = 128. As we see, for high enough temperatures [lines (1), (2) in the figure], 
indeed y{l) is seen to decay very rapidly to zero, LxiP) displaying an exponential behavior. 
However, while for the highest temperature condition [line (1)] the scaling behavior is of 
the /MBE type since Lx(/*) > 1, we observe that there are lower temperatures for which 
the lattice potential is irrelevant, but the long distance behavior is rather of the EW type. 
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since Lx{l*) < 1- This is the case of e.g. the RG flow hne (2) on Figs. ^ and ^ for which 
the system falls onto the high temperature line of flxed points of the sG model. Finally, for 
low enough temperatures the lattice potential dominates the asymptotic properties of the 
system, as signalled in Fig. ^by the complex non-exponential behavior of (/) for T = 0.45. 
These temperature condition corresponds in Fig. ^ to the RG flow of line (3), along which 
the system falls onto the flat low temperature phase of the sG model. In analogy with the sG 
case, the dynamic RG flow (^Bp-(^) thus predicts a roughening transition at a temperature 

T, = M(ii±^^, (44) 

TT 



at which the flow of the lattice potential changes stability. Note both z/ and k in ( ^41) are 
renormalized, rather than bare values. For T < Tr, the system is in the sG massive low 
temperature phase, within which the surface is flat. For T > Tr, the surface is rough and 
should asymptotically feature the EW scale invariant behavior. However, crossover behavior 
exists; speciflcally, for a small system size (for which the asymptotic behavior may not be 
reached) or for very high temperatures [for which the RG flow needs many iterations in order 
to escape the xMBE plane of initial conditions, see line (1) in Fig. 0, the scaling behavior 
observed in the system will be of the /MBE type. 

2. d > 2 cases 

We have also considered the RG flow of the xMBE model for higher substrate dimensions, 
flnding a behavior which is qualitatively consistent with all the above analysis. Moreover, 
given that for any d > 2 the integral (p in ( |i3| ) is always a flnite number, the results that 
follow can be taken as additional support for those presented in ci = 2 employing lattice 
cut-off regularization. As can be seen in Fig. ^, for d > 2 the RG flow of the xMBE model 
does take on the exact shape of the generalized model studied in Sec. [Ill A|, and all the 



conclusions drawn there become applicable also for the speciflc initial condition we are now 
considering. Namely, the RG flow escapes from the x = 2K plane, and thereon the behavior 
is similar to that of the sG model. In particular, there is no phase transition and the surface 
morphology is flat for any value d > 2. Therefore, the upper critical dimension of the 
xMBE model (1)-® is also predicted to be dc = 2. Speciflcally, in Fig. ^, we plot results 
of a numerical integration of (|38D -(^0D starting from the xMBE model initial condition for 
substrate dimensions d = 3 [panel (a)] and d = 4 [panel (b)]. As can be seen in the flgure, 
for any initial condition on the xMBE plane — that is, for any temperature value — , the RG 
flow is towards the sG plane and towards increasing y values, which is the behavior that 
corresponds to the low temperature flat phase, as in the sG model for the corresponding 
substrate dimension. For a given value of d, the higher the temperature is, the longer will it 
take for the RG flow to escape from the plane of initial conditions, since the crossover length 
Lx is larger. In this very high temperature conditions, the initial ZMBE behavior will be 
relevant for a longer time in the dynamics of Eq. (H). 



IV. SUMMARY AND CONCLUSIONS 

The dynamic RG study presented for the generalized model of crystalline surfaces @- 
(0) predicts an upper critical dimension dc = 2. Thus, for a two-dimensional surface there 
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exists a roughening temperature, in such a way that the high temperature phase exhibits the 
scahng properties of the EW equation. Moreover, the transition properties are controlled 
by the sG fixed point (see Appendix and is thus expected to be of a continuous type. If 
we think of (§D-(0) as a generalization of the sG model in which a (possibly small) surface 
diffusion term has been allowed for, we have seen in Sec. Ill A| that the consequence is 
an increase in the sG roughening temperature, due to the crossover induced by such an 
irrelevant perturbation. Moreover, the specific case of the xMBE model has been seen 
to share all these properties, its peculiarity of having a zero bare surface tension merely 
introducing much more severe crossover effects. From the point of view of applications to 
epitaxial growth systems, this result illustrates the relevance of EW scaling as an universality 
class in MBE: in principle, if the symmetry of the system prohibits non-equilibrium surface 
currents and includes invariance of the dynamics under arbitrary surface tilts, ZMBE scaling 
is expected 0]. However, we have obtained for a system relaxing linearly as in the /MBE 
equation that, if a lattice potential influences the dynamics as in (^) — thereby accounting 
e.g. for the discrete character of deposition events or the influence of an underlying lattice — , 
then asymptotic EW scaling should occur. Admittedly, crossovers associated with /MBE 
behavior may be nonetheless rather long, particularly for high temperatures. 

For the case of the xMBE model, the existence of a phase transition between a flat low 
temperature phase and a rough high temperature phase as predicted by (^)-(|40|) is com- 
patible with previous results obtained by the variational approximation [El] and numerical 



simulations JTS], However, while the variational mean fleld predicts the transition to 
be first order, the incorporation of fiuctuations by means of the RG is consistent with a 
phase transition of a continuous type, which is closer to the numerical results. Moreover, 
our present RG study predicts the upper critical dimension to be that of the sG model, 
dc = 2, while mean-field predicts dc = 4. We note that the absence of non-trivial mode 
coupling in the variational mean-field approach limits its capabilities in determining cor- 
rectly both lower and upper critical dimensions, see e.g. |^ and for the case of the 



sG model. Concerning numerical simulations ITR E3|, EW scaling has been obtained in the 



high temperature phase only for temperatures which are extremely close to the roughening 
temperature, whereas /MBE behavior has been obtained for all other temperature values 
above the transition. Recall we obtained in Sec. |1II B 1| that, if the temperature was high 
enough, the system might need a long time to overcome the crossover associated with /MBE 
behavior. Seemingly, the finite system sizes thus far employed in the simulations [L < 128) 
are affected by this type of crossover limitations. We also note that the occurrence of EW 
scaling for the xMBE model is reminiscent of the hexatic phase claimed for the Lr model 



21, 22]. This is an intermediate phase which features EW scaling, and which lies in 



between the flat low temperature phase (the liquid phase in the melting context) and the 
high temperature phase (solid phase) with /MBE scaling, being separated from both of them 



by roughening transitions of the KT type [EH, pTl p2|. On the other hand, as found in ]H5 



and references therein, there are also claims on the inexistence of an hexatic phase in the 
Lr model and, moreover, on the flrst order type of the only roughening transition ensuing. 
We are currently performing large scale simulations of the Lr model on the square lattice 
p9| in order to check the predictions of our RG analysis, in particular whether any trace 
of crossover behavior and EW scaling can be detected in the surface properties at high 
temperatures. 
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APPENDIX A 



In this Appendix we recall the considerations needed to perform the correct Taylor ex- 
pansions in (p!9D, using Kadanoff's operator algebra [0. We generalize results for the static 
RG analysis of the sine-Gordon nonlinearity p|, p3[. Thus, in our dynamical RG calcula- 



tion each of the (infinite number of) non-linearities 02n+i{p) = {h — /i')^"""*"^ appearing in 
the Taylor expansion of sin[27r(/i — h')/a±] contributes a term proportional to the marginal 
operator Oi{p) = h — h', where, as above, p = r — r'. In our continuum approach, this is 
the most relevant term originating (via Taylor expansion) the renormalization of both the 
Ah and A'^h terms in Eq. (M). Thus, 



02„+i(p) = (Mr) - Hr')^'-' 

= 02n+i{r) + a2n+i{p){h - h') 



(Al) 



where 02n+i(r) is an irrelevant operator. A way to compute the a2„+i constants is |32|, |33 
by performing all contractions contributing to the behavior 

(02„+i(r)0i(r + p))h ~ a2„+i(p)(0i(r)0i(r + p))r,. 

where, within our order of approximation in powers of V, averages {■ ■ ■)h are computed with 
respect to the Gaussian distribution of h given by the V = limit of (|^). Thus one obtains 

a2„+i(p) = {2n + mh - h'rh = ^^(Ch " 
Using this result in the Taylor expansion of the sine, one gets 



sm 



2n{h - h') 



I 1 C'2„+l(p) 



(2ra+l)! \a^J 



n=0 

oo / ^ x„ / „ N 2?i+l 



n=0 

27r 
a_L 



2"ra! \a±J 



Oi(r) exp 



271^ 



{{h-h'Yk 



(A2) 



This is the result employed in Eq. ( [20| ) of Sec. ^ 
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APPENDIX B 



As already stated, the dynamic RG flow (|38| ) -(l40| ) can be seen as a generalization of that 
for the sG equation, as derived e.g. in P, In this Appendix we explore this relationship 
in some more detail in the physically interesting case d = 2. 

Similarly to the sG case, an important point in the {x, y, K) parameter space is that 
where the lattice potential y changes stability. Actually, the existence in the generalized 
model of a surface diffusion term does not change this fact, in the sense that the point 
(1, 0, 0) still controls the behavior of the RG flow to a large extent. This is due to the facts 
that the surface diffusion is irrelevant with respect to surface tension, and that the latter is 
generated by the RG flow, as obtained in Sect. The simplest way to substantiate this 
conclusion is to study the flow perturbatively near the point (1, 0, 0). To this end, 

we introduce the temperature-like variable |^ t = 2/a; — 2, rewrite Eqs. (^S|)-(^Dp in the 
new variables (t, K), and approximate the corresponding flow to second order around the 
flxed point {t,y,K) = (0,0,0). We obtain 



dy 

dl 
dK 

'df 



-rr = -ty^ 



-2K 



t)-yHB^'^ 



^r(4) 

64 ' 



(Bl) 
(B2) 
(B3) 



where we have termed i?*-"^ = -^'■"•'(l, 0). Introducing new variables in order to bring flow 
([B1|)-(P3D into normal form 



we flnally obtain 



t 

y 

K 

dui 
IT 

dU2 

'df 

du3 

dl 



Ul + -Ml, 
U2, 



U3 



5(4) 

T28 



-u 



25 



8m3 - 5(2)^2 



-UiU2, 



25 



-2m, 



(B4) 
(B5) 
(B6) 



The interesting feature of the approximate flow (|B4D -( p6D is that the third equation can 
be readily solved for the surface diffusion-like variable M3, as u^i^l) = ^3(0) exp(— 2/), thus 
making apparent one of the qualitative features of the original RG flow, namely the fact 
that surface diffusion is an irrelevant variable that decouples (exponentially) fast from the 
flow for the lattice potential and the surface tension. On the other hand, by taking the ratio 
between the flrst and second equations in ( |B4| )-(p6D and integrating in /, one can check that 
the relation 



ulil) 



/oo 
Uiil')u3il')dl' 



(B7) 
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defines a separatrix for tlie flow ( |B4| )-( p6D , generalizing the well-known asymptotes of the 
sG hyperbolae [see also panel (a) of our Fig. |I|. On e.g. the first quadrant of the (mi,M2) 
plane, flow lines below the separatrix flow onto the U2 = line of fixed points (irrelevant 
lattice potential, high temperature behavior), whereas flow lines above the separatrix flow 
towards large values of the lattice potential U2 (low temperature, massive phase). Thus, the 
separatrix marks the temperature driven phase transition. Unfortunately, we have not been 
able to produce a useful simpler expression for locus (P7|), not even within perturbation 
theory. 
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FIG. 1: Panels (a)-(c): {x,y) projection of the RG flow (|8|)-(||) for d = 2, using u{0) = V{0) = 
1 = a_L = A = 1, and (a) k(0) = 0, (b) k{0) = IQ-"^, (c) k{0) = 10"^ In panels (a)-(c), initial 
conditions lie on the dashed line and the dotted line is the separatrix for the pure sG flow. Solid 
lines correspond to T = 0.5, 0.6, 0.7, 0.8, 0.9, and T = 1 right to left in (a) and (b), and top 
to bottom in (c). Panel (d): Specific heat for the sG model (left curves) and for model @ with 
/i = 1/ = 1 and K = 0.5 (right curves) for different system sizes. All units employed are arbitrary. 
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FIG. 2: RG trajectories from a numerical integration of (p8|)-(^) for d = 3 (a), d = 4 (b), and 
d = 5 (c) with 1/(0) = 0.5, «;(0) = 0.25, and V{0) = 0.1. Thick sohd hnes correspond, right to 
left, to (a) T = 1.5, 2, 2.5, 3, 3.5, 10; (b) T = 10, 15, 20, 30, 40; (c) T = 5, 10, 15, 20, 25. Thin 
solid lines on the (x, y) planes are projections of the RG flow lines above them. Thick dots on the 
X axis denote both the origin and the point x = 2SciA'^~'^ /[d{27rY^^] at which the rhs of the flow 
equation for y vanishes. Other parameters are as in Fig. [l|. All units employed are arbitrary. 
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FIG. 3: RG trajectories from a numerical integration of (^8|)-(^) for the xMBE model [i.e. ^{O) = 0] 
with L = 128, k(0) = 1, and V{0) = 0.1. Solid lines exemplify the three types of behavior depending 
on T: (1) high T phase, /MBE scale invariant behavior (T = 10); (2) intermediate T phase, EW 
scale invariant behavior (T = 0.46); (3) low T massive phase (T = 0.45). For the sake of clarity, 
all three coordinates along line (1) have been artificially expanded by a factor of 10. The x = 2K 
plane of initial conditions appears shaded, and both the origin and the (1, 0, 0) point are signalled 
with thick dots. Other parameters are as in Fig. |l|. All units employed are arbitrary. 




FIG. 4: Linear-log plot of the numerical RG flow of Lx{l) for d = 2, L = 128. Other parameters 
are as in Fig. ^. Solid lines are for T = 10, 0.46, and T = 0.45, top to bottom. For the sake of 
comparison, the dashed line depicts the RG flow of Lx{l) for the linear system V{0) = 0. In all 
cases the flow is terminated for I = In 128 ~ 4.85. All units employed are arbitrary. Inset: Blow-up 
of the same plot for small I values, in log-log representation. 
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FIG. 5: RG trajectories from a numerical integration of (p8|)-(pO|) for the xMBE model [i.e. z^(0) = 0] 
with k(0) = 1, and V{0) = 0.1. Panel (a) [(b)] corresponds to d = 3 [d = 4]. On each panel, solid 
lines correspond to temperature values T = 1.5, 2, 2.5, 3, 3.5, and T = 10, right to left. Thin 
solid lines on the (x, y) plane are projections of the RG flow lines above them. As in Fig. |^, the 
x(0) = 2i^(0) plane of initial conditions appears shaded. Thick dots on the x axis denote both the 
origin and the point x = 2Sd^'^''^ /[d{2Tr)'^^'^] at which the rhs of the flow equation for y vanishes. 
Other parameters are as in Fig. ||. All units employed are arbitrary. 
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